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Beschreibung 
Technischer Hintergrund 

5 Im Bereich der chemischen Industrie gelten die sogenannten Life Science Produkte als der erfolgversprechendste Be- 
reich der nachsten Jahre. Da jedoch pharmazeutische Praparate und Feinchemikalien immer kornplexeren Anforderun- 
gen und strengeren gesetzlichen Auflagen gerecht werden mussen, miissen hier effiziente Verfahren zur schonenden Pro- 
dukttrennung, wie z. B. chromatographische Trennverfahren, zum Einsatz kommen. Der kontinuierliche Betrieb chro- 
matographischer Trennverfahren mit scheinbarem Feststoffgegenstrom als SMB-ProzeB (Simulated-Moving-Bed, Si- 
10 mulierter Gegenstrom) gewinnt zur Trennung einer StofFmischung in zwei Fraktionen zunehmend an Bedeutung. Wegen 
der hohen Produktkosten stellt der wirtschaftliche Betrieb dieser Anlagen hohe Anforderungen an die ProzeBfuhrung 
UDd ProzeBiiberwachung. 

Stand derTechnik 

15 

Der SMB-ProzeB wird durch die Hintereinanderschaltung einer beliebigen Anzahl einzelner Russigchromatographie- 
saulen oder -segmente erreicht, und zwar so, daB sich eine endlose Schleife ergibt. Ublich sind Anordnungen mit 8 bis 24 
Saulen oder Segmenten. Edukt- und Losungsmittelzufuhr sowie Abfuhr der beiden Produkte werden zwischen den Sau- 
len angeordnet und definieren vier verschiedene Zonen des Prozesses, die jeweils zwischen den Zu- und Ablaufen liegen. 

20 Der scheinbare Gegenstrom wird naherungsweise durch ein periodisches Umschalten der Zu- und Ablaufe in Richtung 
des Ruidstroms realisiert, d. h. nach einer spezifizierten Schaltperiode werden alle Zu- und Ablaufe jeweils eine Saule 
weitergeschaltet. Zone I liegt zwischen der Desorbenszufuhr und der Extraktabfuhr, Zone H zwischen der Exraktabfuhr 
und der Feedzufuhr, Zone HI zwischen Feedzufuhr und Raffinatabfuhr und Zone IV zwischen Raffinatabfuhr und Desor- 
benszufuhr. Die Zuordnung der Zonen erfolgt aufgrund ihrer Funktionalitat und ihrer relativen Lage zu den Zu- und Ab- 

25 laufen. Die Zonen werden nicht bestimmten Chromatographiesaulen zugeordnet, die Funktion der einzelnen Saulen an- 
dert sich mit jeder Schaltperiode. Eine schematische Skizze des Prozesses befindet sich in Abb. 1. Adsorptive Trennver- 
fahren nach diesem Prinzip sind beschrieben in (Broughton et al. 1961), (Rasche, 1992), (Negawa and Shojii, 1995), 
(Priegnitz, 1995). Detaillierte Beschreibungen der moglichen technischen Ausfuhrungen dieses Prozesses konnen den 
voranstehendenReferenzen entnommen werden. Sie dienen auBerdem derDokumentation des Standes der Technik. Eine 

30 technische Ausfuhrung des Prozesses in einer einzelnen Saule ist beschrieben in (Schoenrock, 1992), in (Hotier, 1998) 
wird eine Betriebsweise mit veranderlichen Schaltperioden und konstantem Kreislaufstrom beschrieben. Der SMB-Pro- 
zeB zeichnet sich durch eine komplexe, aus regelungstechnischer Sicht schwer zu beherrschende, gemischt kontinuier- 
lich/diskrete Systemdynamik aus. Problematisch sind dabei vor allem der Anfahrvorgang und der stabile Betrieb der An- 
lagen in der Nahe des wirtschaftlichen Optimums. Die derzeit bekannten Verfahren zur Regelung von SMB-Prozessen 

35 beruhen uberwiegend auf einer reinen Vorgabe der relevanten Russigkeitsstrbme. In (Masao and Tanimura, 1986) wird 
eine Methode beschrieben, diese Russigkeitsstrome innerhalb der Saulen zu regeln. Dieses Verfahren zielt darauf hin, 
die Riissigkeitsstrome mit Hilfe von Druckmessungen an einem konstanten Sollwert zu halten. Ein ahnliches Verfahren 
ist in (Cohen et al., 1997) beschrieben. Ein Verfahren mit einem iiber die Zeit variablen, aber fest vorgegebenem, Soll- 
wert beschreibt (Kearney and Hieb, 1992). Diese Verfahren beruhen allesamt nicht auf der Messung der tatsachlichen 

40 Systemzustande und Produktkonzentrationen, sondem benutzen ausschlieBlich die im voraus berechneten Volumen- 
strome, so daB aufgrund der in der Umgebung des optimalen Betriebspunktes hohen Sensitivitat des Prozesses gegeniiber 
Storungen und Parameterschwankungen ein stabiler Betrieb der Anlagen in der Nahe des wirtschaftlichen Optimums 
nicht immer moglich ist. 

Neueste Untersuchungen haben gezeigt, daB fur die Regelung des SMB-Prozesses nicht allein die Produktkonzentra- 

45 tionen wesentlich sind, sondem daB auch die Lage und die zeitliche Entwicklung samtlicher Adsorptions- bzw. Desorp- 
tionsfronten fur die ProzeBfuhrung von entscheidender Bedeutung sind. Diese GroBen sind jedoch nicht direkt meBbar, 
da entsprechende Konzentrationsmessungen innerhalb der Saulen technisch nicht moglich sind und auch MeBstellen 
nach jeder einzelnen Saule oft einen unvertretbaren finanziellen Aufwand bedeuten. Erschwerend kommt hinzu, daB die 
ProzeBdynamik durch variable, von den Betriebsparametern abhangige Totzeiten bestimmt wird und daB - je nach Posi- 

50 tion der Konzentrationsplateaus in den Saulen - die Sensitivitat der QualitatsgroBen beziiglich der EinfluBgroBsen sehr 
stark variiert. Dies fuhrt unter anderem auch dazu, daB Storungen oft erst nach sehr langer Zeit (mehrere Stunden bis 
Tage) sichtbare Auswirkungen haben. Bis heute gibt es nur wenige Ansatze, die konventionelle Fahrweise durch eine 
fortgeschrittene Regelung unter Einbeziehung von Konzentrationsmessungen undBeobachtung der intemen Konzentra- 
tionsprofile zu verbessem, um so die Prozesse in der Nahe des wirtschaftlichen Optimums zu betreiben. In (Holt, 1995a) 

55 ist ein Regelungs- und MeBkonzept fur einen SMB Prozess beschrieben, welches in (Holt, 1995b) noch einmal speziell 
fur eine para-xylene Trennung entwickelt wurde. In diesem Verfahren werden aus on-line Messungen der Konzentratio- 
nen an zwei Stellen des Kreisiaufstromes mit Hilfe eines Algorithmus die Riissigkeitsstrome und Umschaltperioden so 
bestimmt, daB eine Produkteigenschaft, wie z. B. die Reinheit, einen Sollwert einhalt. Im Gegensatz zu der hier beschrie- 
benen Erfindung erfordert dieses Verfahren weitere MeBstellen zur Produktiiberwachung und erlaubt auch keine Beob- 

60 achtung der Konzentrationsprofile in den Saulen. 

In (Marteau et al., 1994) ist ein Kpnzept zur ProzeBiiberwachung einer Aromatentrennung mittels SMB Chromatogra- 
phic beschrieben. Die dort vorgeschlagene MeBstellenanordnung enthalt vier MeBstellen, wobei zwei im Kreislaufstrom 
angeordnet sind und jeweils eine die Konzeritration des Extrakt- und Raffinatstroms iiberwacht Dieses Verfahren beno- 
tigt also vier MeBstellen fur eine "Oberwachung des Prozesses und der Produktzusammensetzung. Aus Kostengriinden 

65 ware es jedoch vorteilhaft, weniger MeBstellen zu verwenden. 

Die hier beschriebene Erfindung soli die Nachteile der obigen Ansatze vermeiden. Insbesondere soil sie ermoglichen, 
mit nur zwei MeBstellen sowohl die Produktqualitat on-line zu uberwachen als auch den Verlauf des Konzentrationspro- 
fils innerhalb einer Saule verfugbar zu machen. Aus der Kennmis des Konzentrationsprofils in den Saulen lassen sich 
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dann vor allem die zur ProzeBfuhrung wesentlichen Positionen der Adsorptions- und Desorptionsfronten ableiten. Auf 
dieser Grundlage kann z, B. im weiteren ein Regeialgorithmus entwickelt werden, mit dem die geforderten Produktei- 
genschaften auch im Falle von Storungen und Parameterschwankungen eingehalten werden. 

Die hier dargestellte Erfindung bildet zum einen dafiir die Grundlage, stellt aber gleichzeitig auch eine Methode dar, 
mit der das nichtmefibare Konzentrationsprofil der Anlage zu jeder Zeit uberwacht werden kann sowie die relevanten 5 
ProzeBparameter on-line ermittelt werden konnen, 

Beschreibung der Erfindung 

Die Erfindung betrifft eine Methode zur Parameterschatzung und ProzeBbeobachtung fur simulierte Gegenstrom- 10 
Chromatographieprozesse. Die Erfindung erlaubt die kontinuierliche Ermittlung der in alien Sauien vorliegenden Kon- 
zentrationsprofile aus der kontinuierlichen Messung der AusfluBkonzentration an zwei Stellen des Prozesses, 

Die Erfindung ist besonders geeignet zur ProzeBiiberwachung von simulierten Gegenstrom-Chromatographieprozes- 
sen im Bereich der Zuckertrennung (Rearick et al., 1997) z. B. Aufspaltung von Fructose-Glucose Losungen, Abtren- 
nung von Sucrose aus Molasse, (Broughton, 1983), 15 

Die MeBeinrichtungen werden direkt in die Leitungen fiir die beiden Produkte Extrakt und Raffinat installiert und lie- 
fern on-line die Konzentrationen der beiden Komponenten, Eine fur eine Zuckertrennung taugliche MeBapparatur ist 
z. B. beschrieben in (Altenhoner et al., 1997). 

Bei der Anordnung der MeBstellen muBte beriicksichtigt werden, daB die Eigenschaften der Chromatographiesaulen 
von der Geometrie und Packung der Sauien abhangen, welche aus produktionstechnischen Griinden leichten Schwan- 20 
kungen unterliegen. Im Extremfall hat jede Einzelsaule also voneinander abweichende Eigenschaften. Die MeBstellen 
wurden deshalb so angeordnet, daB aus den MeBdaten die Systemparameter saulenspezifisch ermitteln konnen, Aus Ko- 
stengriinden ist es ein weiteres Anliegen, moglichst wenige MeBstellen zu verwenden. Die hier gewahlte Anordnung der 
MeBstellen hat zusatzlich den Vorteil, daB die im Rahmen der Profil- und Parameterermittlung durchgefuhrten Messun- 
gen gleichzeitig zur Produktuberwachung genutzt werden konnen. 25 

Der Erfindung liegt zugrunde, daB fur eine Vielzahl von Trennungen das physikalische Verhalten der Adsorptionssaule 
hinreichend genau durch folgende Gleichung beschrieben werden kann: 

dci dd d 2 Ci 

Hierbei ist q(x,t) die Konzentration der Komponente i am Ort x zur Zeit t und u die Stromungsgeschwindigkeit in den 
Sauien. Eine Beschreibung dieses Modells, seiner physikalischen Grundlagen und und seiner Ldsung findet sich in (Dim- 
nebier et al., 1998). Yi und D^ (i = A, B) sind vier fur das System charakteristische Parameter, welche aus den MeBwer- 
ten on-line bestimmt werden miissen. 35 

Parameter- und Profilermittlung 

Unter der Voraussetzung, daB die Substanzen A und B bei der Adsorption keine Wechselwirkungen aufeinander aus- 
iiben, konnen die Parameterpaare y A und D^ sowie Yb und D^ voneinander getrennt ermittelt werden. 40 

Die nachfolgend beschriebene Parameterermittlung ist jeweils saulenspezifisch. Diejenige Saule, deren Parameter be- 
stimmt werden, wird als "relevante Saule" bezeichnet. 

Die in Extrakt- bzw. Raffinatstrom gemessene AusfluBkonzentration c- umcss (x = L,t) bei unbekannten Parametem yi 
und Dap^ (i = A, B) wird mit einer numerisch ennittelten AusfluBkonzentration q(x = L, t; % T>apJ) abgeglichen. In der 
numerischen Simulation sind die Parameter yi und Dap^ bekannt Der Abgleich erfolgt fiber diese beiden Parameter durch 45 
Optimierung des quadratischen Giitefunktionals 

J(7i,D* P ,i)= (2) 

J (c i , me ^(t)-Ci(x = L,t;ji,D a . p ,i)^ dt so 

(i = A,B) . 

iiber einen noch zu spezifizierenden Zeithorizont (ti, ti). 55 
q(x = L, t; % D^O ist die Losung der partiellen Differentialgleichung (1) mit Rand- und Anfangsbedingungen 



ci(x = 0,t; Yi,Dap4) = Cc4(t) (3) 
Ci(x,t = 0; yi, Dap.0 = C 0 ,l(x). (4) 



60 



Mit den bzgl. (2) optimalen Parametem Yi = Yi,opt und ^>aps = Dap4,opt gibt Ci (x,t = 0; Yuopt, Dap j,opt) das Konzentrations- 
profil am Taktbeginn im Inneren der relevanten Saule wieder. 

Die Rand- und Anfangsbedingungen (3) und (4) sind nicht bekannt und werden approximiert. Zur Ermittlung der 
Rand- und Anfangsbedingungen werden die gemessenen Konzentrationen in Extrakt- und Rafiinatabzug in zwei aufein- 
anderfolgenden Takten benotigt. Die Parameter konnen nach dem zweiten Takt bestimmt werden. Dieser Takt wird des- 65 
halb als "aktueller Takt" bezeichnet Zur Unterscheidung wird der erste Takt, aus dem MeBungen verwendet werden, als 
"vorangegangener Takt" angesprochen. 

L bezeichnet die Saulenlange und te die Taktzeit Als "Richtung" wird die FlieBrichtung des Desorbens gewahlt 
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Bestininiung von y A undD^ 

Die Betriebsparameter der Anlage miissen so gewahlt werden, dafi (nach einer Anlaufzeit) die Desorptionsfront der 
Substanz A, die die schwacher adsorbierbare Komponente bezeichnet, den Raffinatentnahmepunkt EX in jedem Takt 
5 passiert. Nach der Anlaufpbase bildet sich im Bereich des Extraktabzugs in zwei aufeinanderfolgenden Takten die in 
Abb. 2 gezeigte Konzentrationsverteilung aus. 

Zur Bestirnmung der Parameter y A und D^ und des Konzentrationsprofils werden die Auslaufkonzentrationen im 
Extraktabzug des vorangegangenen und des aktuellen Taktes, c^me^aitW und CAjnessW, verwendet. 

10 Approximation der Anfangsbedingung 



Im vorangegangenen Takt ist die relevante Saule die erste.Saule von Abschnitt II (vgl. Abb. 2a). Die Einlaufkonzen- 
tration in diese Saule entspricht der Konzentration im Extraktabzug und ist somit mefibar. 

In der zweiten Zone bildet sich im Konzentrationsprofil der Substanz A ein Plateau annahemd konstanter Konzentra- 
15 tion aus (Diinnebier et al„ 1998). Aus diesem Grund wird das Anfangsprofil als naherungsweise konstant gleich 
C A ,ni e ss,ak(t = 0) gewahlt 

Zur Approximation des Konzentrationsprofils in der relevanten Saule am Ende des vorangegangenen Taktes, c Kal £x t t 
= te; y A , D^), wird die partieUe DGL (1) mit Rand- und Anfangswerten 



20 C efAta it( t ) c AfinesSiaIt (t) 
gelost. 

Nach dem Umtakten ist die relevante Saule die letzte Saule von Abschnitt I (vgl. Abb. 2b). Somit entspricht der An- 
25 fangswert im aktuellen Takt 

co^W = c A^t(x,t = te;Y A ,D aftA ) (5) 

30 Approximation der Randbedingung 

Aus grober Kenntnis des ungefahren Profils ist klar, daB die Konzentration der Substanz A in Zone I zu Ende eines je- 
den Taktes sehr gering ist. Nennenswerte Konzentrationen finden sich nur im Endbereich der letzten Saule. Nach dem 
Umtakten ist die vormals letzte Saule von Abschnitt I die vorletzte Saule von Abschnitt I (vgl. Abb. 2b und c). Aus die- 
35 sem Grund wird die nicht meBbare Einlaufkonzentration in die relevante Saule im aktuellen Takt im wesentlichen als 
Null angenommen. Der Anfangsbereich wird durch die empirisch ermitteite Formel 

CeAt)= ( 6 ) 

n eS s,ait(te)(l - tl.85j A ) *e [0,1.857a] 
0 sonst 



40 | CA,n 



linear approximert. y A istdabei der aktuelle Wert der Optimierung. (6) entspricht einer Modellierung des Konzentrations- 
45 profils durch lineare Approximation und anschlieBender Approximation des Ausflusses durch ein Modell ohne Diffu- 
sion. 

Die Losung von (1) mit Rand- und Anfangswerten (6), (5) approximiert c A (x,t; Ya, Vapud in der relevanten Saule. 
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Festlegen des Zeithorizonts 

Der Zeithorizont des.Gutefunktionals (2) wird so gewahlt, daB die Desorptionsfront der Substanz A den Extraktabzug 
innerhalb dieses Zeitraums passiert 

Bestirnmung von Yb und D^ 

Die Betriebsparameter der Anlage miissen so gewahlt werden, dafi (nach einer Anlaufzeit) die Adsorptionsfront der 
Substanz B den Raffinatentnahmepunkt RA in jedem Takt passiert (vgl. Abb. 3). 

Zur Bestirnmung. der Parameter Yb "nd D^ und des Konzentrationsprofils werden die Auslaufkonzentrationen im 
Extraktabzug des vorangegangenen und des aktuellen Taktes, c^mwitfl) und CB,mess(0> verwendet. 

Approximation der Anfangsbedingung 



Die Anfangsbedingung Co,b(x) wird unter Verwendung der Einlaufkonzentration der Substanz B in die relevante Saule 
im vorangegangenen Takt, CB tm cas^it(0» approximiert Im vorangegangenen Takt ist die relevante Saule die erste Saule 
65 von Abschnitt IV. Die Einlaufkonzentration in diese Saule ist mefibar, da sich der Raffinatabzug im vorangegangenen 
Takt vor der betrachteten Saule befindet (vgL Abb. 3a). 

CB,mess.alt(0 = CB^nessW (7) 
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Aus grober Kenntnis des ungefahren Profils ist klar, daB die Konzentration der Substanz B in Zone IV zu Beginn eines 
jeden Taktes sehr gering ist (vgl. Abb. 3a und c). Nennenswerte Konzentxationen finden sicb nur im Anf angsbereich der 
jeweils ersten Saule. Aus diesem Grund wird die Losung der partiellen Differentialgleichung (1) in der relevanten Saule 
(vorangegangener Takt) approximiert, indem als Randverteilung zum Zeitpunkt t = 0 im wesentlichen die Nullinie ange- 
nommen wird. Der Anfangsbereich wird durch die empirisch ermittelte Formel 

Cg,bM x ) ( 8 ) 



f c B , mesS|alt (0)(l-t5/L) *€[0,L/5] 
| 0 sonst 



Zusamraenfassung 

Als Ergebnis liefert das Verfahren mit Hilfe von nur zwei MeBstellen on-line: 
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Patentanspriiche 

20 1. Eine Methode zur on-line Parameterschatzung und ProzeBbeobachtung fur SMB (simulierte Gegenstrom-) 

Chromatographieprozesse, dadurch gekennzeichnet, daB mit Hilfe von zwei kontinuierlichen Messungen, jeweils 
an den Produktentnahmestellen angeordnet, aus den MeBdaten mit Hilfe eihes mathematischen Algorithmus die das 
physikalische Verhalten des Systems bestimmenden Parameter in Echtzeit bestimmt werden. 
2. Ein Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB aus der MeBinformation und mit Hilfe einer Be- 

25 rechnungsvorschrift die innerhalb der Chromatographiesaulen vorliegenden nicht meBbaren Konzentrationen ge- 
schatzt werden. 
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Abbildung 1: Prinzipskizze des SMB-Prozesses 
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Abbildung 2: Ungefahrer Verlauf des nicht mefcbaren Konzentrationsprofils von Substanz A 
im Bereich der Desorptionsfront wahrend zweier aufeinanderfolgender Takte 
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Abbildung 3: Ungefahrer Verlauf des nicht mefibaren Konzentrationsprofils von Substanz B 
im Bereich der Adsorptionsfront wahrend zweier aufeinanderfolgender Takte 
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